Topological defects: fossils of an anisotropic era? 
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We consider the evolution of domain walls produced during 
an anisotropic phase in the very early universe, showing that 
the resulting network can be very anisotropic. If the domain 
walls are produced during an inflationary era, the network will 
soon freeze out in comoving coordinates retaining the imprints 
of the anisotropic regime, even though inflation makes the 
universe isotropic. Only at late times, when the typical size 
of the major axis of the domain walls becomes smaller than 
the Hubble radius, does the network evolve rapidly towards 
isotropy. 

Hence, we may hope to see imprints of the anisotropic era 
if by today the typical size of the major axis of the domain 
walls is of the order of the Hubble radius, or if the walls re- 
entered it only very recently. Depending on the mass scale 
of the domain walls, there is also the possibility that they 
re-enter at earlier times, but their evolution remains friction- 
dominated until recently, in which case the signatures of the 
anisotropic era will be much better preserved. These effects 
are expected to occur in generic topological defect models. 

PACS number(s): 98.80. Cq, 95.30.St 



I. INTRODUCTION 

It is well know that the 'Hot Big Bang' model Q, de- 
spite its numerous successes, is plagued by a number of 
'initial conditions' problems, of which the horizon, flat- 
ness and unwanted relic ones are the best known. The 
standard way to solve them is to invoke an epoch of cos- 
mological inflation j^i} , a relatively brief period of ex- 
ponential (or quasi-exponential) cosmological expansion. 
The way inflation solves these problems is, loosely speak- 
ing, by erasing all traces of earlier epochs and re-setting 
the universe to a rather simple state. Indeed, inflation 
is so efficient in this task that a number of people have 
wondered if one can ever hope to probe the physics of a 
pre-inflationary epoch. 

There are, however, a small number of possible pre- 
inflationary relics. For example, the recent work by 
Turok and collaborators [^) shows that curvature can, 
in some sense, survive inflation. Another class of infla- 
tionary survivors are topological defects ||, formed at 
phase transitions either before or during inflation 
It is known (see, e.g. l|||[lCfl|) that one needs about 20 
e-foldings of inflatiorj^pto solve the monopole problem. 
One can reverse the argument and say that monopoles 
can survive about 20 e-foldings of inflation. The infla- 
tionary epoch itself will obviously push the monopoles 
outside the horizon, but the subsequent evolution of the 
universe tends to make them come back inside, so if the 
inflationary epoch is not too long they can still have im- 
portant cosmological consequences. 

Cosmic strings are even more successful, being able to 
survive about 50 e-foldings. The reason for this differ- 
ence is that their non-trivial dynamics [[fl] P"3||ic| ] makes 
them come back inside the horizon faster than one might 
naively have expected. The above two numbers are typ- 
ical, but there are specific models where defects can sur- 
vive even longer. One example is that of open inflation 
scenarios In this case the universe undergoes two 

different inflationary epochs — roughly speaking, a period 
of 'old inflation' followed by one of 'new inflation'. As 
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pointed out by Vilenkin p5[ , one can expect that defects 
will form between the two inflationary epochs. In this 
case, a collaboration jlO| including the present authors 
has recently shown that not only will cosmic strings sur- 
vive the entire second inflationary epoch, regardless of 
how long it lasts[| but they will in fact be back inside 
the horizon by the time of equal matter and radiation 
densities. In such models, monopoles can survive up to 
about 30 e-foldings. 

Now, given that defects seem to be so successful sur- 
viving inflation, and that one expects them to be frozen 
out while they are outside the horizon, one can think of a 
further interesting possibility. For the best-studied case 
of cosmic strings, it is well known that the scaling proper- 
ties of the network depend on the background cosmology 
p| p^[ . Moreover, in some cases (typically when their 
evolution is friction-dominated) they can retain a 'mem- 
ory' of the initial conditions, or the general properties 
of the cosmology in which they find themselves at early 
times, for quite a large number of orders of magnitude 
in time 0,[13|]. It is therefore conceivable that if such 
an imprint of an early cosmological epoch is retained by 
a defect network which manages to survive inflation, we 
might still be able to observe it today. 

We believe that this is a general feature of defect mod- 
els, and a number of non-trivial pieces of information 
about the very early universe can probably be preserved 
in this way. In the present paper we will restrict our- 
selves to a simple example. We will discuss the possibil- 
ity of a domain wall network retaining information about 
an early anisotropic phase of the universe. There are 
very strong constraints [fl6|| on the mass of domain walls 
formed after inflation, due to the fact that their density 
decays more slowly than the radiation and matter den- 
sities. However, these can be evaded by walls forming 
before or during inflation. In a subsequent paper, we will 
discuss the more interesting, but also more complicated, 
case of cosmic strings. 

The plan of the paper is as follows. In section || we 
briefly describe our background (Bianchi I) cosmology 
and the basic evolutionary properties of the domain walls. 
In particular, we focus on the approach to isotropy dur- 
ing inflation, which is discussed through both analytic 
arguments and numerical simulations. We emphasise 
that these simulations do not include the defects. How- 
ever, they serve an important purpose, as they are used 
in the subsequent discussion to show that the timescale 
needed for isotropization is compatible with the 'survival' 
on anisotropic defect networks. 

We provide a description of our numerical simulations 
of domain wall evolution in section III. These arc analo- 



2 Note that in these models the duration of the second infla- 
tionary epoch is fixed by the present value of the density of 
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gous to those of Press, Ryden & Spergel [0, and the 
interested reader is referred to this paper for a more 
detailed discussion of some relevant numerical issues. 
Here defect networks are evolved in an isotropic, matter- 
dominated (ie, post-inflationary universe), and their 
main purpose is to show that isotropic and anisotropic 
networks will evolve in different ways, so two such net- 
works can in principle be observationally distinguished as 
they re-enter the horizon. Our main results are presented 
and discussed in section IV, and finally we present our 
conclusions and discuss future work in section 

Throughout this paper we will use fundamental units 
in which c = I. 



II. EVOLUTION EQUATIONS FOR DOMAIN 
WALLS 

We consider the evolution of a network of domain walls 
in a k=0 anisotropic universe of Bianchi type I with line 
element [H: 



ds 2 = dt 2 - X 2 (t)dx 2 - Y 2 (t)dy 2 - Z 2 (t)dz 2 



(1) 



where X(t),Y(t) and Z(t) are the cosmological expansion 
factors in the x, y and z directions respectively, and t is 
the physical time. The dynamics of a scalar filed cf> is 
determined by the Lagrangian density, 



where we will take V{4>) to be the generic 
with two degenerate minima given by 



v (<l>) = M 72 - 1 



(2) 

i> 4 potential 
(3) 

§. Byvary- 
(4) 



This obviously admits domain wall solutions 
ing the action 



S = J dt J d^xyf^gC, 
with respect to </> we obtain the field equation of motion: 

(5) 
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with 0(f) = W/W and W(t) = XYZ. The dynamics of 
the universe is described by the Einstein field equations. 
Here we shall seek perfect fluid solutions. The time com- 
ponent of the Einstein equation then becomes 



+ A 2 + B 2 + C 2 = -\k{p + 3p), 



(7) 
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while the spatial components give 



A + 0A = B + 6B = C + 6C= -k(p-p), (8) 

with A = X/X, B = Y/Y and C = Z/Z, 9 = A + B + C 
and k = 8ttG/c 2 and i = 1, 2, 3. It is straightforward to 
combine equations to obtain: 

(9) 



AB + BC + CA — kp. 



In the following discussion we will make the simplification 
that X(t) — Z(t) (and therefore A = C) and consider the 
dynamics of the universe during an inflationary phase 
with p = —p = const. In this case H 2 = kp/3 = const 
and the Einstein field equations imply: 



A+^(A 2 - H 2 ) = 0, 
while B can be found from the suggestive relation 



B 
~A 
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A 2 
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(10) 



(11) 



Equation (|T^) has two solutions, depending on the initial 
conditions. If A; < H, then A is the smaller of the 
two dimensions and the shape of spatial hyper-surfaces 
is similar to that of a rugby ball. Then the solution is 



A 

— = tanh 
H 



^H(t-U 



tanh" 1 | — 



(12) 



with A; = A{ti). On the other hand, if Aj > H, then A 
is the larger of the dimensions and the shape of spatial 
hyper-surfaces is similar to that of a pumpkin. Then the 
solution is 



A 
H 



= coth 



-H(t - U) + coth" 



(13) 



Note that in both cases the ratio A/H tends to unity 
exponentially fast, and hence the same happens with the 
ratio B I A. In other words, inflation tends to make the 
universe more isotropic, as expected. An easy way to see 
this is to consider the ratio of the two different dimen- 
sions, D = B/A, and to study its evolution equation. 
One easily finds 



1/2 



(1--D), 



(14) 



which has an obvious attractor at D = 1. 

Note that even though we have so far assumed (for 
simplicity) that p = —p, the same analysis can be carried 
out for an inflating universe with p = (7— l)p with 7 7^ 
by numerically solving the conservation equation 



p + 9(p + p) = 0, 



(15) 



together with equations (||) and (|^). Indeed, the more 
general case will be relevant for what follows. 



In figure [l] we plot the evolution of the asymmetry pa- 
rameter E — Y(t)/X(t), according to eqns. @j and (|), 
for several values of on — log 10 (Di) assuming 7 = and 
7 = 2/3 (note that 7 = 2/3 is the maximum value of 7 
which violates the strong energy condition). Note that 
we do not include the defect network in the simulation. 
(We assume that the network at the initial time U is sta- 
tistically isotropic.) We take X(ti) — Y(ti). We can see 
that depending on the initial degree of anisotropy, spec- 
ified by on, the value of E can grow to be very large, 
especially if on is large. Moreover, although for 7 = 0, 
the value of E becomes approximately constant in one 
Hubble time that does not happen so rapidly for inflat- 
ing universes with larger 7. This removes the necessity 
of producing the domain walls right at the onset of the 
inflationary era. 

What about the evolution of the domain walls? Based 
on rather general grounds, we expect it to have a number 
of similarities with the much better studied case of cosmic 
strings ||. In particular, one can define a 'characteristic 
length-scale', which we shall denote by L, that can be 
roughly interpreted as a typical curvature radius or a 
correlation length of the wall network. It is also a length- 
scale that measures the total energy of the domain wall 
network per unit volume, since we can define 



P =L 



(16) 



where a is the domain wall energy per unit area. Note 
that in a more rigorous treatment that allowed for the ex- 
pected build-up of small-scale 'wiggles' on the walls (in 
analogy with what happens for the case of cosmic strings 
p3|) each of these three length scales would be differ- 
ent. However, for our present purposes it is adequate to 
suppose that they are all similar. 

Then we can expect to find two different evolution 
regimes. While the network is non-relativistic, we ex- 
pect it to be conformally stretched by the cosmological 
expansion, and hence 



L cc a . 



Pw eta 



(17) 



In this case there is essentially no dynamics. An extreme 
example of this regime happens during inflation We can 
see from eqn. (^|) that due to the very rapid expansion 
which occurs in the inflationary regime the time deriva- 
tives of the field (j> rapidly approach zero so that the net- 
work of domain walls will simply be frozen in comoving 
coordinates. 

On the other hand, once the network becomes rela- 
tivistic, one expects it to evolve in a linear scaling regime 
where 



L oc t . 



Pw OC t 



(18) 



This is the case of 'maximal' dynamics, in the sense that 
the network is evolving (in particular, losing energy by 
wall collisions and re-connections) as fast as allowed by 
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FIG. 1. Evolution of the asymmetry parameter E — Y(t)/X(t) (according to eqns. (Q) and (@)) for several values of 
Oi = log 10 (Di) assuming 7 = (left panel) and 7 = 2/3 (right panel) respectively. Note that time is given in units of if (tj). 



causality. We note that previous work of Press, Ryder 
and Spergel [[l7| suggests that there may be logarithmic 
corrections to this linear regime. 



III. NUMERICAL SIMULATIONS 

At late times (after the inflationary epoch) the universe 
is homogeneous and isotropic with A = B = C with 
the average dynamics of the universe being specified by 
the evolution of the scale- factor a(t). We now consider 
the evolution of isotropic and anisotropic defect networks 
in this background. In particular, we are interested in 
determining how the networks evolve as they re-enter the 
horizon, since if one finds differences in the dynamics of 
the two cases then this should translate into observational 
tests that will allow us to discriminate between then and 
hence probe pre-inflationary physics. 

It is useful for numerical purposes to re- write equation 
(^) as a function of the conformal time r\ defined by drj = 
dt/a. In this case equation (|J) becomes 



d 2 <j) ^adcj) 
drj 2 a drj 



X7 2 d> = 



,dV 



with 



dx 2 



dy 2 



dz 2 



(19) 



(20) 



When making numerical simulations of the evolution of 
domain wall networks (or indeed other defects) it is also 
often convenient to modify the equation of motion for the 
scalar field dj in such a way that the comoving thickness of 
the walls is fixed in comoving coordinates. This is known 
as the PRS algorithm |L7j , and it is generally believed not 
to significantly affect the large-scale dynamics of domain 
walls. 

We note, however, that recent high-resolution simula- 
tions |HJ have revealed that the accuracy of this algo- 
rithm is not as good as has been claimed. This effect 
is expected to increase with increasing dynamic range. 
In particular, the PRS algorithm artificially prevents the 
build-up of small-scale features on the domain walls (or, 
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for that matter, any other defect). This turns out to be 
crucial for a quantitatively accurate description of their 
evolution, and hence for a reliable analysis of their obser- 
vational consequences. For our purposes in the present 
work, however, the PRS algorithm is enough as an ap- 
proximation to the true wall dynamics. In a subsequent, 
more detailed publication we shall compare results ob- 
tained using this algorithm with those from the true wall 
dynamics. 

Having clarified this point, we will modify the evolu- 
tion equation for the scalar field <f> in the isotropic phase 
according to the PRS prescription: 



a i 

a drj 



where (3i and are constants. We choose 02 = in 
order for the walls to have constant comoving thickness 
and 0i = 3 by requiring that the momentum conservation 
law for how a wall slows down in an expanding universe 
is maintained [ |j"7| |. 

We perform two-dimensional simulations of domain 
wall evolution for which d 2 cf>/dz 2 — 0. These have the 
advantage of allowing a larger dynamic range and better 
resolution than tree-dimensional simulations. 

We solve equation (|2l]) numerically assuming a matter- 
dominated Einstein-de-Sitter cosmology with a cx rf . We 
used a standard difference scheme second-order accurate 
in space and time and periodic boundary conditions (see 
p7[ for a more detailed description of the algorithm and 
other related numerical issues). 

The initial properties of the network of domain walls 
depend strongly on the details of the phase transition 
which originated them. It is conceivable that the initial 
network is already formed asymmetric with the walls be- 
ing elongated along preferred directions. However, this is 
beyond the scope of the present paper. For our present 
purposes, we can ignore this possibility and assume that 
the initial domain wall network is statistically isotropic. 
This assumption will not modify the conclusions of the 
paper — if anything, any ab initio anisotropics would only 
enhance the effects we are describing. 

Hence, we assume the initial value of <j) to be a random 
variable between — </>o and 4>o and the initial value of 4> to 
be equal to zero everywhere. We normalise the numerical 
simulations so that 4>o — 1- We set the conformal time 
at the start of the simulation and the comoving spacing 
between the mesh points to be respectively r\i = 1 and 
Ax = 1. 

The wall thickness, defined by 



UJ 



7T0Q 



(22) 



is set to be equal to 5. The kinetic energy of the field 
is calculated by 



On the other hand, the rest energy of the walls is calcu- 
lated by multiplying the comoving area of the walls, A, 
by the energy density per comoving area, which can be 
written as 



2~f w V uj 
3tt 2 ' 



(24) 



(23) 



with 7m = (1— v^,) -1 ' 2 , and v w is the value of the physical 
velocity of the domain walls. Finally the total area of the 
walls is defined as the area of the surfaces on which </> = 
and is computed using the method described in ref . |l7j . 

IV. RESULTS AND DISCUSSION 

As pointed out above, it will be of fundamental im- 
portance to study the dynamics of the wall network at 
late times, as the Hubble length becomes larger than 
the typical size of the major axis of a domain wall. A 
crucial issue will be the timescale required for the wall 
network to switch from the non-relativistic regime to 
the relativistic one. For our present purposes, the main 
difference between these two regimes is that a friction- 
dominated network can remain anisotropic if it froze out 
that way, whereas a relativistic network will rapidly be- 
come isotropic and erase any imprints from the earlier 
anisotropic phase. In our simulations we ignore the pos- 
sibility that the network can be friction dominated due 
to particle scattering Jl3| when the domain walls come 
back inside the horizon — again, this would only enhance 
the effects we are describing. 

We consider three simulations with different initial 
conditions. In the first one (case I) we evolve the initial 
network generated in the manner specified in the previous 
section from the conformal time rji — 1. In the second 
one (case II) the initial conditions at the time rji = 1 
were specified by the network configuration of the previ- 
ous simulation at the conformal time 77* =20, with the 
velocities reset to zero. Physically, this corresponds to 
starting with the network outside the horizon. Finally, 
the case III is similar to the second one but with the 
initial network of case II stretched in the y direction by 
a factor of E = 2 (see fig. ||), and corresponds to the 
anisotropic case. We have performed 1024 2 simulations 
for each of the three cases, plus an additional 2048 2 run 
of case I, in order to test for possible box effects. 

For each run we plot (see fig. ||) the ratios A/V and 
Aq/V (note that A and V are the comoving area and 
volume, respectively), as well as the ratio of the kinetic 
and rest energies, as in Press, Ryden and Spergel [[L7| . 
These are plotted from the beginning of the simulation 
until the time when the horizon becomes one half (for 
the 1024 2 runs) or one quarter (for the 2048 2 run) of the 
box size. In addition to these (which we plot mainly for 
the purposes of comparison with previous work p7|) we 
plot a 'scaling coefficient' which will be our main analysis 
tool. We will define it by analogy with the cosmic string 
case [0,[l3[, as follows. Assume that 
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FIG. 2. The same physical size of the same domain wall simulation. The bottom one however has been stretched along the 
y direction by a factor of 2. The horizon size of the top box is 1/8 of the linear size of the box shown, and this is itself only a 
fraction 1/16 of the side of the whole simulation box.] 
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(25) 



then what we plot is the 'instantaneous' or 'effective' 
value of A as a function of conformal time. For a given A, 
the physical network correlation L length will be evolving 
as 



L 



L(xt 1 ' X/3 , -(XT] 1 



— A 



(26) 



Note that A can, in general, be a time-dependent quan- 
tity. However, for the two scaling regimes discussed 
above, we expect it to be a constant, namely 



Arir 1 



(27) 



in the non-relativistic limit where the network is being 
conformally stretched, and 



A r = 



(28) 



in the linear scaling regime. 

From these it is trivial to deduce the behaviour of A/V 
and Arj/V in both scaling regime. One expects 



in agreement with the previous result by Press, Ryden 
and Spergel |L7|]. This means that the network is not 
straightening out as fast as allowed by causality. Sec- 
ondly, the rates at which the networks in cases II and 
III approach the relativistic regime are different. One 
might expect this on physical grounds: if the network is 
stretched in one direction, then there are in fact differ- 
ent 'network correlation lengths' for each direction, and 
interactions between the domain walls will tend to oc- 
cur faster along the directions with smaller correlation 
lengths, and more slowly in the others. 

Another way of saying this is that the network will 
only start evolving towards the relativistic regime when 
its larger axis has re-entered the horizon. Note that this 
mechanism also tends to make the domain wall network 
more isotropic. So one can naively say that the approach 
to the linear scaling regime takes longer in an anisotropic 
universe because the dynamics of the walls must accom- 
plish two tasks (make the wall network relativistic and 
isotropic) rather than just one. 



V. CONCLUSION 
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(29) 



(30) 



in the linear scaling regime. Similarly, the ratio E^/ E rest 
should be a constant in the linear scaling regime (with its 
numerical value providing a measure of the characteristic 
network scaling speed), and it should approach zero in 
the non-relativistic limit. 

Firstly, we note that the two case I runs produce very 
similar results: significant differences can only be seen 
at late times. This is an indication that the resolution 
we are using is adequate for our present purposes. As 
expected, the network in case I becomes relativistic very 
quickly, while those of cases II and III start in the extreme 
non-relativistic regime and only evolve away from it fairly 
slowly, after they re-enter the horizon. 

More importantly, there are two non-trivial observa- 
tions to be made. Firstly, we confirm that there is a 
correction to the linear scaling regime. We find 



A sc ~0.12, 



(31) 



which corresponds to to evolve in a linear scaling regime 
where 



L oc t 



l0.96 



Pw oc i 



-0.96 



(32) 



In this paper we have discussed a simple example of 
what we believe to be a rather generic feature of topolog- 
ical defect models, namely that they can easily retain in- 
formation about the properties of the very early universe. 
This information is encoded in the scaling (ie, 'macro- 
scopic') and statistical (ie, 'microscopic') properties of 
the defect networks. This is even more relevant given 
the fact that defects can survive significant amounts of 
inflation. Hence, they can provide a unique probe of the 
pre-inflationary universe. The two crucial scales in the 
problem are the defect mass scale and the epoch when 
the defects come back inside the horizon. 

Specifically, we have discussed the role of domain walls. 
We have highlighted the existence of two scaling regimes 
for the domain wall network, in agreement with previ- 
ous work jl7|]. Furthermore, we have shown that an 
anisotropic network re-entering the horizon will take 
longer to approach scaling than an isotropic one. Hence, 
if the very early universe had an anisotropic phase which 
was erased by an inflationary epoch, and if domain walls 
are present, then the walls can retain an imprint of the 
earlier phase, and this can have important observational 
consequences, eg for structure formation scenarios. 

As is well known, there are quite strong constraints 
jl6|,H on the mass of domain walls formed after inflation. 
These are basically due to the fact that their density will 
decay more slowly than the radiation and matter densi- 
ties. However, essentially all of these can be evaded (or 
at least significantly relaxed) by walls forming before or 
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during inflation (and also by walls evolving in a friction- 
dominated regime). Having said this, how could these 
anisotropics be detected? The most naive answer would 
be through their imprint on CMB, but this is only true 
if their energy density is not too low, and such models 
are constrained in a variety of other ways (not only from 
the cosmology side, but also from the high-energy physics 
side). The case of 'light' walls is therefore more interest- 
ing: note that just like in the case of 'light strings' [ p^[ , 
these are expected to be friction-dominated throughout 
most of the cosmic history. Here the observational de- 
tection of the effects we have described becomes some- 
what non-trivial. The best way of doing it should be 
through observations of numbers of objects as a func- 
tion of redshift in different directions (assuming that one 
has a reliable understanding of other possible evolution- 
ary effects). Two specific examples would be large-scale 
velocity flows |2(J and gravitational lensing statistics of 
extragalactic surveys |21| . 

Finally, there is also an important implication of our 
work if at least one of the minima of the scalar field poten- 
tial has a non-zero energy density, which is an anisotropic 
non-zero vacuum density. In a subsequent, more detailed 
publication, we shall discuss this scenario in more detail, 
as well as the analogous one for cosmic strings. 

To conclude, we have shown that the importance of 
topological defects as a probe of cosmological physics 
goes well beyond structure formation. Even if defects 
turn out to be unimportant for structure formation they 
can still (if detected) provide us with extremely valu- 
able information about the physical conditions of the very 
early universe. 
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FIG. 3. Evolution of several properties of a domain wall network as a function of the conformal time rj during the matter 
era. Here, Ek and E re st are the kinetic and rest energies associated with the scalar filed <f>, A/V is the comoving wall area per 
unit comoving volume of the two-dimensional simulations, A is a 'scaling coefficient' defined by A = d\n(Ari /V) / dint] and r] is 
the conformal time. The solid, dashed and dash-dotted curves correspond to 1024 2 simulations with different initial conditions. 
Cases I and II have isotropic initial conditions but in case two the size of the domain walls becomes comparable with horizon 
only at a conformal time r\* = 20. Case III is similar to case II but with the initial network of case II stretched in the y direction 
by a factor of E — 2 (see text). The dotted curve corresponds to a 2048 2 run of case I. 
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